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ABSTRACT 


,  CWf« /Ey-P^-olf  fi-) 

During  the  ICE^flyby  of  Glacoblnl-Zlnner  the  plasma  wave  instru¬ 
ment  detected  numerous  impulsive  signals  caused  by  dust  Impacts  on  the 
spacecraft.  Most  of  the  Impacts  occurred  within  30,000  km  of  the 
comet.  The  Impact  rate  averaged  over  the  Inbound  and  outbound  legs 

r  SfC'/ir-c, 

varies  approximately  as  l/;t^ ,  as  would  be  expected  for  an  isotropic 
constant  velocity  radial  outflow.  Small  differences  between  the 
Inbound  and  outbound  legs  exist  which  may  be  indicative  of  azimuthal 
variations  in  the  dust  production  rate.  A  simple  model  of  the  Impact 
ionization  and  charge  collection  by  the  antenna  Indicates  that  the  par¬ 
ticles  have  a  mass  on  the  order  of  to  10“^ °  gm,  corresponding  to 

particles  with  radii  of  a  few  microns.  '  •  '  '  -  / 
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I.  INTRODUCTION 


During  the  ICE  (International  Cometary  Explorer)  flyby  of  the 
comet  Glacoblnl-Zlnner  on  September  11,  1985,  the  plasma  wave  Instru¬ 
ment  detected  numerous  Impulsive  electric  field  signals  that  are 
believed  to  be  due  to  dust  Impacts  on  the  spacecraft  [Scarf  et  al., 
1986].  The  purpose  of  this  paper  Is  to  present  an  Initial  analysis  of 
these  Impacts. 

Dust  Impact  effects  of  this  type  were  first  observed  during  the 
Voyager  2  crossing  of  Saturn's  ring  plane.  During  the  ring  plane 
crossing  both  the  plasma  wave  and  radio  astronomy  instruments  detected 
numerous  Impulsive  noise  bursts  on  the  electric  antenna  [Scarf  et  al., 
1982;  Warwick  et  al.,  1982).  Subsequent  analysis  by  Gurnett  et  al. 
[1983]  and  by  Aubler  et  al.  [1983]  led  to  the  conclusion  that  the 
Impulsive  signals  were  caused  by  mlcron-slzed  particles  striking  the 
spacecraft.  When  a  small  particle  strikes  a  solid  surface  at  a  suffi¬ 
ciently  high  velocity,  the  particle  is  Instantly  vaporized  and  Ionized. 
Laboratory  measurements  show  that  the  charge  released  Is  directly 
proportional  to  the  mass  of  the  Impacting  particle  [for  a  review,  see 
Fechtig  et  al.,  1978].  Some  of  the  charge  released  by  the  Impact  is 
collected  by  the  electric  antenna,  thereby  producing  a  voltage  pulse. 
The  amplitude  of  the  voltage  pulse  is  therefore  proportional  to  the 
mass  of  the  particle. 


II.  OVERVIEW  OF  OBSERVATIONS 


An  expanded  16-channel  plot  of  the  antenna  voltages  detected  dur¬ 
ing  a  representative  2-minute  interval  near  the  closest  approach  to 
Giacobini-Zinner  is  shown  in  Fig.  1.  The  output  of  each  channel  is 
proportional  to  the  logarithm  of  the  voltage  amplitude  with  a  dynamic 
range  of  90  db,  from  about  10  pVolt  to  0.3  Volt.  The  output  of  each 
channel  is  processed  by  a  peak  detector  which  discharges  exponentially 
with  a  time  constant  of  0.3  sec.  The  16  channels  are  sampled  sequen¬ 
tially  starting  with  the  100-kHz  channel  and  ending  with  the  17,8-Hz 
channel.  The  time  required  for  a  scan  of  all  16  channels  is  0.125  sec, 
and  a  complete  spectrum  scan  is  obtained  once  per  second.  For  a 
further  description  of  the  Instrument,  see  Scarf  et  al.,  [1978]. 

As  can  be  seen  from  Fig.  1,  the  antenna  voltages  are  characterized 
by  numerous  impulsive  events,  often  extending  simultaneously  across  a 
broad  range  of  frequencies.  Some  of  these  events  have  a  well-defined 
peak  in  the  spectrum,  typically  at  frequencies  from  about  1  to  50  kHz. 
The  events  at  1107:23,  1107:39  and  1107:49  are  of  this  type.  We  inter¬ 
pret  these  broadband  bursts  as  plasma  waves.  As  described  by  Scarf  et 
al.  [1986],  numerous  impulsive  plasma  wave  emissions  of  this  type  were 
observed  near  the  comet. 
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Another  type  of  Impulsive  event,  with  a  more  nearly  monotonic 
spectrum,  can  be  seen  in  Fig.  1  decreasing  in  intensity  with  Increasing 
frequency.  The  events  at  1106:31,  1106:43  and  1107:15  are  of  this 
type.  We  Interpret  these  events  as  dust  impacts.  Spectrums  of  these 
events  are  shown  in  Fig.  2.  These  events  are  believed  to  be  dust 
impacts  because  they  are  qualitatively  similar  to  the  Impacts  detected 
by  Voyager  2  at  Saturn.  The  power  law  index,  y,  of  the  spectrum  varies 
from  about  -2  to  -4,  with  a  tendency  for  the  spectrum  to  steepen  at 
higher  frequencies,  very  similar  to  the  Voyager  2  spectrum  given  by 
Aubler  et  al.  [1983].  The  ICE  spectrums  are  not  as  Intense  as  the 
Voyager  2  spectrums  because  the  impact  rate  is  not  as  high  (the 
spectrums  in  Fig.  2  are  single  Impacts,  whereas  the  spectrum  given  by 
Aubier  et  al.  is  an  average  over  many  Impacts),  but  otherwise  are  quite 
similar. 


V  V 


in.  IMPACT  RATE  AND  NUMBER  DENSITY 


In  order  to  determine  the  impact  rate,  a  computer  algorithm  was 
developed  to  identify  impulsive  events.  Because  of  the  large  dynamic 
range,  the  identification  criterion  was  based  on  the  logarithm  of  the 
voltage  spectral  density,  I  =  logjQ(V^/Af).  Impulsive  events  were 
selected  by  requiring  that  the  difference  between  preceding  and  follow¬ 
ing  intensities,  AI_  =  1^  “  In-1  ^1+  “  In  “In+1»  averaged  over  the 
56.2  Hz  through  10-kHz  channels,  exceed  a  preset  threshold.  Channels 
below  56.2  Hz  were  not  used  because  of  the  high  Interference  levels  in 
these  channels,  and  channels  above  10  kHz  were  not  used  because  the 
impulses  are  usually  below  the  Instrument  noise  levels  at  these  fre¬ 
quencies.  A  threshold  value  of  AIj^^  =  0.2  was  found  to  give  a  reliable 
identification  of  impulsive  events  without  introducing  a  large  number 
of  spurious  events. 

To  separate  impacts  from  plasma  wave  emissions  a  power  law  fit 
test  was  utilized.  Again  the  56.2-Hz  through  10-kHz  channels  were 
used.  Events  were  classified  as  impacts  if  the  power  law  index  (slope 
of  the  best  straight  line  fit  of  1  versus  logjg  f)  was  in  the  range 
-4.0  <  Y  <  -2.5,  and  the  r.m.s.  error  of  the  fit  was  less  than  0.7. 

This  power  law  fit  test  was  found  to  give  very  good  discrimination 
against  plasma  waves  because  the  spectrum  of  plasma  waves  typically  has 


a  peak,  which  gives  a  very  poor  fit  to  a  power  law.  Using  this  impact 

identification  criterion  a  total  of  294  Impacts  were  identified  during 

the  2-hour  period  from  1000  to  1200  UT. 

Next  we  consider  the  counting  rate  as  a  function  of  time.  For 
this  analysis  it  is  useful  to  introduce  an  amplitude  threshold,  so  that 

only  particles  with  masses  greater  than  a  certain  threshold  are  count¬ 

ed.  The  r.m.s.  amplitude  of  the  voltage  pulse  can  be  computed  from  the 
area  under  the  spectrum,  Vj-^g^  =  /(V^/Af)df.  The  counting  rate  R  is 
shown  in  Fig.  3  for  two  thresholds,  1  mVolt,  and  10  mVolt.  The  error 
bars  give  the  statistical  uncertainty  in  the  counting  rate.  At  closest 
approach  the  counting  rates  range  from  about  2  to  10  counts  per  minute 
depending  on  the  threshold.  The  rates  for  the  highest  threshold,  10 
mVolt,  are  probably  the  most  reliable,  because  they  represent  the 
largest  signals.  However,  even  for  1  mVolt,  which  is  near  the  back¬ 
ground  noise  level,  the  count  rate  shows  a  well-defined  peak  centered 
on  closest  approach,  indicating  that  these  events  are  probably 
reliable . 

The  counting  rate  has  also  been  computed  as  a  function  of  radial 
distance  from  the  comet,  and  is  shown  in  Fig.  4.  Because  of  the  low 
counting  rate  only  impacts  greater  than  1  mVolt  are  considered,  and  the 
inbound  and  outbound  legs  have  been  averaged  together  to  Improve  the 
counting  statistics.  The  results  of  this  analysis  show  a  smooth 
monotonic  decrease  with  increasing  radial  distance.  The  slope  of  the 
counting  rate  profile  on  the  log-log  plot  is  roughly  consistent  with  a 
1/r^  density  variation,  which  is  what  would  be  expected  for  an 
Isotropic  radial  outflow  of  dust  from  the  comet  at  a  constant  velocity. 


The  number  density,  n,  can  be  computed  from  the  counting  rate 
using  the  relation 

R  =  nUA  (x/T)  ,  (1) 

where  U  is  the  relative  speed  between  the  particles  and  the  spacecraft, 
A  is  the  detecting  area,  and  x/T  is  a  dead  time  correction  due  to  the 
discharge  of  the  peak  detector.  Since  the  velocity  of  the  dust  rela¬ 
tive  to  the  comet  is  small  compared  to  the  spacecraft  velocity,  U  is 
simply  the  approach  velocity  which  is  about  21  km/ sec.  The  dead  time 
correction  can  be  computed  from  the  discharge  time  constant  of  the  peak 
detector,  x  =  0.3  sec,  and  the  time  between  samples,  T  =  1  sec,  which 
gives  x/T  =  0.3.  The  largest  uncertainty  in  the  number  density  arises 
from  the  area  factor.  The  detecting  area  depends  on  whether  the  rele¬ 
vant  impacts  are  on  the  spacecraft  body  or  the  antenna.  The  flyby  tra¬ 
jectory  is  such  that  the  particles  arrive  at  an  angle  of  about  12.8® 
from  the  +Z  axis  of  the  spacecraft,  roughly  parallel  to  the  spin  axis. 
Most  of  the  impacts  are  therefore  on  the  top  of  the  spacecraft,  which 
has  an  area  of  about  2.4  m^ .  For  Voyager  2  it  has  been  shown  [Gurnett 
et  al.,  1983]  that  the  relevant  impacts  were  on  the  spacecraft  body. 
Therefore,  based  on  this  consideration  alone,  the  best  choice  for  the 
detecting  area  would  be  A5(,  =  2.4  m^.  A  scale  giving  the  number  den¬ 
sity  computed  using  A  =  2.4  m^  is  shown  on  the  right  side  of  Fig.  4. 

Two  factors  make  us  believe  that  the  detecting  area  may  be  less 
than  2.4  m^  ,  First,  using  Ag^,  =  2.4  m^  the  peak  flux  and  total  number 
of  impacts  are  several  orders  of  magnitude  below  the  pre-encounter 


values  given  by  Yeomans  and  Brandt  [1985],  based  on  a  dust  model  by 
Divine  et  al.  [1985],  Second,  inspection  of  the  spacecraft  shows  that 
the  top  is  recessed  well  below  the  sides  (like  a  can  with  the  lid 
removed),  so  that  the  antennas  are  shielded  from  Impacts  on  the  top  of 
the  spacecraft,  thereby  reducing  the  effective  area.  The  minimum 
possible  detecting  area  is  the  area  of  the  antenna.  Since  impacts  on 
the  antenna  deposit  charge  directly  on  the  antenna,  all  antenna  impacts 
should  be  detected.  For  particles  arriving  from  the  +Z  direction,  the 
projected  area  of  the  antenna  is  =  0.053  m^ .  The  number  density 

using  this  area  is  shown  by  the  scale  marked  A  =  0.053  m^  on  the  right 
side  of  Fig,  4.  These  number  densities  are  larger  and  in  better  agree¬ 
ment  with  the  pre-encounter  predictions  of  Yeomans  and  Brandt,  Since 
Asc  3nd  A^j^t  bound  the  possible  detecting  areas,  the  two  scales  on  the 
right  of  Fig.  4  give  upper  and  lower  limits  to  the  number  density,  n^j^^ 


IV.  MASS  DISTRIBUTION 


Next  we  consider  the  difficult  problem  of  estimating  the  mass  of 
the  particles  from  the  amplitude  of  the  voltage  pulse.  From  laboratory 
measurements  it  is  known  that  the  charge  emitted  Is  proportional  to  the 
mass  m  of  the  particle, 

Q  =  km,  (2) 

where  k  is  a  yield  constant  that  depends  on  the  speed  and  composition 
of  the  particle.  For  an  impact  speed  of  21  km/sec  typical  yield  con¬ 
stants  for  dielectric  particles  like  silicon  dioxide  or  water  ice  are 
about  k  =  1  Coul/gm  [Dalmann  et  al.,  1977].  Variations  of  up  to  a 
factor  of  10  can  occur  for  various  particle  compositions  and  target 
materials.  We  will  use  a  yield  constant  of  1  Coul/gm.  The  amplitude 
of  the  voltage  pulse  on  the  antenna  can  then  be  written 

V.«(iS),  (3) 

where  is  the  capacity  of  th-?  antenna  (including  base  capacity)  and  a 
is  a  coefficient  that  takes  into  account  the  fraction  of  the  emitted 
electrons  collected  by  the  antenna.  For  the  45m  antenna  element  on  ICE 
the  capacity  C^  is  about  265  pf  (personal  communication,  R.  Stone, 
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Because  of  the  complicated  geometry  of  the  spacecraft  and  the 
unknown  bias  voltage  of  the  antenna,  it  is  difficult  to  accurately 
estimate  the  collection  coefficient  a.  The  maximum  value  is  obviously 
“max  “  following  calculations  we  will  adopt  a  =  1.  Equation 

3  can  then  be  used  to  place  a  lower  limit  on  the  particle  mass, 


m  >  m*  = 


In  fact,  m*  is  probably  close  to  the  actual  mass  because  estimates  of 
the  collection  coefficient  for  Impacts  on  Voyager  2  [Gurnett  et  al., 
1983J  ,  give  collection  coefficients  close  to  one  (a  =  0.6).  For 
impacts  on  the  antenna,  the  collection  coefficient  should  be  exactly 


One  final  correction  must  be  made.  Because  the  pulse  duration  is 


short  compared  to  the  response  time  of  the  filters,  the  pulse  ampli¬ 
tude,  V,  is  actually  larger  than  the  r.m.s.  voltage,  Vj-^g,  computed 
from  the  spectrum.  For  a  short  pulse  the  correction  factor  can  be 


shown  to  be 


where  Tq  is  the  duration  of  the  pulse,  and  At  =  l/nAf  is  the  filter 


response  time.  From  the  Voyager  waveform  measurements  a  typical  pulse 
duration  is  known  to  be  a  few  tenths  of  one  msec  [Gurnett  et  al.. 


1983].  Since  most  of  the  contribution  to  Vj-^g  occurs  near  100  Hz  (see 
Fig.  2),  At  is  taken  to  be  the  response  time  at  100  Hz.  For  100  Hz  the 


Cl? 


filter  bandwidth  is  Af  =  15  Hz,  so  At  =  21  msec.  Assuming  a  typical 
pulse  duration  of  Tq  *  0.5  msec,  the  correction  factor  for  the  filter 
response  is  approximately  V  =  6.5  Vj.jj,g.  Although  the  correction  factor 
depends  somewhat  on  the  duration  and  shape  of  the  pulse,  the  error 
Involved  is  probably  small  compared  to  the  uncertainties  in  k  and  a. 

Based  on  the  above  values  the  threshold  voltages  of  1  mVolt  and  10 
mVolts  used  in  the  previous  section  are  found  to  correspond  to  masses 
of  m*  =  1.6  X  10“^^  and  1.6  x  10~^^  gm.  Assuming  a  mass  density  of  1 
gm/cm^ ,  these  masses  correspond  to  particles  with  radii  of  0.7  and  1.6 
pm.  Given  that  the  actual  mass  is  probably  somewhat  larger  than  the 
computed  minimum  mass ,  we  conclude  that  the  Impacts  are  produced  by 
particles  with  a  radii  of  a  few  microns. 

Because  the  particle  mass  is  proportional  to  the  pulse  amplitude, 
the  pulse  amplitudes  provide  a  direct  measurement  of  the  mass  distribu¬ 
tion.  The  pulse  amplitude  distribution,  dN/dVj.^s,  integrated  over  the 
Interval  from  1000  to  1200  UT,  is  shown  in  Fig.  5.  A  scale  giving  the 
mass,  m*,  computed  using  the  previously  mentioned  parameters,  is  shown 
at  the  top  of  the  diagram.  Corresponding  particle  radii  are  also 
given,  assuming  a  density  of  1  gm/cm^.  This  plot  shows  that  the  mass 
distribution  has  a  peak  between  about  10~^^  and  10”^^  gm,  and  then 
decreases  rapidly  with  increasing  mass,  varying  approximately  as  dn/dm 
~  m“^ .  The  amplitude  distribution  near  1  mVolt  is  sensitive  to  the 
choice  of  the  counting  threshold  Al^^*  so  the  exact  shape  of  the  dis¬ 
tribution  in  this  region  (indicated  by  a  dashed  line)  is  somewhat 


uncertain 
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V.  DISCUSSION 


A  strong  case  has  been  presented  that  the  impulsive  antenna 
voltages  detected  during  the  ICE  flyby  of  Giacobini-Zlnner  are  caused 
by  dust  impacts  on  the  spacecraft.  The  radial  distribution  of  the 
impact  rate  is  roughly  consistent  with  the  expected  1/r^  density  varia¬ 
tion  for  an  isotropic  radial  outflow  at  a  constant  velocity,  and  the 
size  of  the  particles  determined  from  the  pulse  amplitudes  is  con¬ 
sistent  with  infrared  evidence  for  the  existence  of  micron-sized  dust 
grains  near  a  comet  (Manner,  1984;  Divine  et  al.,  1985],  The  greatest 
uncertainty  in  the  present  analysis  is  in  the  detecting  area.  With 
further  study  we  hope  to  reduce  the  uncertainty  in  some  of  the  param¬ 
eters  Involved,  possibly  by  combining  these  measurements  with  suitable 
optical  and  Infrared  data. 

Although  the  radial  variation  of  the  Impact  rate  is  generally  con¬ 
sistent  with  a  radial  outflow  model,  significant  differences  do  exist 
between  the  inbound  and  outbound  legs.  For  example,  the  abrupt  step  in 
the  count  rate  at  1045  UT  on  the  inbound  leg  does  not  appear  at  a 
comparable  radial  distance  on  the  outbound  leg.  These  differences  are 
probably  due  to  azimuthal  asymmetries  in  the  dust  production  rate. 
Radiation  pressure  also  should  affect  the  dust  distribution  near  the 
comet.  Analysis  of  these  and  othe..  effects  will  require  further  study. 
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FIGURE  CAPTIONS 

Fig.  1.  A  16-channel  plot  of  the  electric  antenna  voltages  near 
closest  approach.  The  sharp  spikes  extending  across  all 
channels  from  about  31  Hz  to  10  kHz  are  believed  to  be  dust 
impacts. 

Fig.  2.  The  spectrum  of  some  selected  dust  Impacts  from  Fig.  1. 

Fig.  3.  The  counting  rate  for  dust  impacts  exceeding  two  thresholds,  1 
and  10  mVolt.  Most  of  the  Impacts  occur  within  30,000  km  of 
closest  approach. 

Fig.  4.  The  radial  distribution  of  the  counting  rate  and  number 
density  averaged  over  the  Inbound  and  outbound  legs. 

Fig.  3.  The  differential  pulse  amplitude  distribution.  This 

distribution  is  proportional  to  the  mass  distribution,  dn/dm. 
A  scale  giving  the  approximate  mass  and  particle  radius  is 


shown  at  the  top  assuming  a  collection  coefficient  of  one 
(a  =  1). 
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